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Synchrotron emission is commonly found in relativistic jets from active galac- 
tic nuclei (AGNs) and microquasars, but so far its presence in jets from young 
stellar objects (YSOs) has not been proved. Here, we present evidence of polar- 
ized synchrotron emission arising from the jet of a YSO. The apparent magnetic 
field, with strength of ~0.2 milligauss, is parallel to the jet axis, and the polar- 
ization degree increases towards the jet edges, as expected for a confining helical 
magnetic field configuration. These characteristics are similar to those found in 
AGN jets, hinting at a common origin of all astrophysical jets. 



Supersonic jets are observed to emerge from a wide variety of astrophysical systems, 
from young brown dwarfs to active galatic nuclei (AGNs). Despite their different physical 
scales (from hundreds to billions of astronomical units), they have strong morphological 
similarites (1), and a common feature of these systems is the presence of a gaseous disk 
around the central object. However, it is yet unclear whether there is a universal mechanism 
that can explain the origin of all these jets. 

It is believed that the formation of stars takes place with the assistance of an accretion 
disk that transports gas and dust from the envelope of the system to the protostar and a jet 
that removes angular momentum and magnetic flux, allowing the accretion to proceed (2). 
Theoretical models for the formation of jets from young stellar objects (YSOs) assume that 
the initial acceleration takes place relatively close to the protostar as gas from the accretion 
disk is lifted and centrifugally accelerated along magnetic field lines. This mechanism was 
first suggested as the origin of jets from accretion disks around black holes (3,4) and was 
soon proposed as the mechanism for protostellar jets (5,6). Present observational techniques 
cannot resolve and directly probe into the innermost regions where the jet is launched, but 
theoretical models for protostellar jets are able to explain some bulk observational properties, 
such as the correlation between accretion and ejection rates (7). The collimation of the 
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outflow into a jet takes place at larger distances from the star and is also believed from 
theoretical arguments to take place magnetically. The magnetic field lines anchored in the 
rotating accretion disk are proposed to wrap in a helical configuration that is needed to 
provide the magnetic pressure to collimate the jet via "hoop" stresses (8,9). Observations 
aimed at measuring the magnetic field in protostellar jets are very scarce and imply complex 
modeling and/or a large number of assumptions (10,11). 

A very powerful method to simultaneously obtain the structure and strength of the 
magnetic field in jets is through observations of the synchrotron emission. When electrons 
move at relativistic velocities in the presence of a magnetic field, they emit linearly polarized 
synchrotron radiation. The spectrum of the synchrotron emission in the centimeter wave- 
length range is characterized by a negative spectral index (i.e., the flux density decreases 
with frequency), and it is related to the strength of the magnetic field. Additionally, the 
direction of the polarization is related to the direction of the magnetic field. This method 
has been the most successful in the study of the magnetic field in jets from AGNs (12) and 
microquasars (13), where the ejected material moves at relativistic velocities. In contrast, 
protostellar jets, where the bulk of the material moves at much smaller velocities (200 to 1000 
km s _1 ) (14,15), apparently would lack relativistic electrons that can provide information 
on the magnetic field. Then, the emission is usually dominated by a different mechanism 
(free-free emission originated from the thermal motions of the electrons, characterized by a 
positive spectral index and no polarization). Intringuingly, in a few protostellar jets (asso- 
ciated with both low- and high-mass YSOs) radio emission with negative spectral indices, 
which could correspond to nonthermal synchrotron emission from relativistic electrons, has 
been reported [Serpens (16), HH 80-81 (17), Cepheus A (18), W3(H 2 0) (19), L778 VLA 
5 (20), IRAS 16547—4247 (21)] (22). If relativistic electrons were present, the acceleration 
of particles would be most likely produced where the fast thermal jet impacts on the sur- 
rounding medium and a strong shock wave is formed. Thus, particles might be accelerated 
up to relativistic energies by the Fermi mechanism (23). In order to confirm the presence of 
synchrotron radiation in these protostellar jets, a clear detection of linearly polarized radio 
emission is required. Because polarized emission is only a fraction of the continuum emission 
and protostellar jets are usually very weak, highly sensitive observations are required. 

We obtained highly sensitive (total observation time ~12 hours) radio continuum images 
of the HH 80-81 jet using the Very Large Array (VLA) in its C configuration (24). This 
highly collimated radio jet is driven by the protostellar source IRAS 18162—2048, whose 
total luminosity is ~17,000 L Q (where L is the luminosity of the Sun) [if this luminosity 
originated in a single object, it would correspond to a zero-age main-sequence star of ~10 M 
(where M & is the mass of the Sun), according to the stellar evolutionary models of (25)]. The 
radio jet consists of a chain of radio sources aligned in a southwest-northeast direction. To 
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the southwest, it terminates in the optical and radio Herbig-Haro objects HH 80 and HH 81, 
whereas to the northeast it terminates in the radio Herbig-Haro object HH 80 N (17) (Fig. 
1). With a total extension of 5.3 pc [for an assumed distance of 1.7 kpc (26)], this is one of 
the largest and most collimated protostellar jets known so far. Previous radio observations 
showed that the spectrum of the emission from the central source is characterized by a 
positive spectral index, suggesting that it is dominated by thermal free-free emission (17). 
In contrast, the negative spectral indices of the emission from HH 80-81, HH 80 N, as well 
as from some of the knots in the jet lobes, suggest that an additional nonthermal component 
could be present in these sources (17). 

Our observations (Figure 1) show that the emission of the knots located ~0.5 pc from 
the driving source is linearly polarized, indicating the presence of nonthermal synchrotron 
emission in this jet (24) and implying the presence of relativistic electrons and a magnetic 
field. The observed polarization vectors are perpendicular to the direction of the jet, with a 
degree of polarization of the order of 10 to 30%. The direction of the apparent magnetic field 
(the component in the plane of the sky averaged over the line of sight) is taken to be per- 
pendicular to the polarization vectors [this is correct for a nonrelativistically moving source, 
such as the HH 80-81 jet; for relativistic jets, such as AGN jets, additional assumptions on 
the velocity field are required (27)}. Then, the apparent magnetic field appears very well 
aligned with the direction of the HH 80-81 radio jet (Fig. 1C). 

By using the equations 7.14 and 7.15 of (28) we obtained the minimum total en- 
ergy, E = Cis(l + A;) 4 / 7 (f) 3 / 7 R 9 / 7 Lr 7 , and the equipartition magnetic field, B = [4.5ci2(l + 
In these equations, <p is the filling factor of the emitting region, R is the 
radius of the source, L R is the integrated radio luminosity, k is the ratio between the heavy 
particle energy and the electron energy, and Ci 2 and Ci 3 are functions of the spectral index 
and of the minimum and maximum frequencies considered in the integration of the spectrum, 
which are given in appendix 2 of (28). To estimate the radio luminosity, we integrated the 
radio spectrum in the range between 20 and 2 cm by using the flux density measurements of 
(17) (table SI). We used a filling factor of 0=0.5 and A;=40 [this value of k is appropriate 
for acceleration in a nonrelativistic shock (29)}. With these parameters, we obtained typical 
values for the knots of the radio jets of B ~ 0.2 mG and E ~ 2 x 10 43 erg. 

The direction of the apparent magnetic field obtained from the synchrotron emission is 
strongly concentrated in a direction very close (~ 4°) to the jet axis, whereas the direction of 
the magnetic field inferred from polarization of the dust near the protostar (30) is offset by 
~ 20° and shows more scatter (Figure 2); however, part of this scatter could be attributed 
to turbulent and thermal motions. It has been estimated that the magnetic field in the 
region ~0.1 pc around the protostar IRAS 18162-2048 has a value of 0.2 mG (30). If we 
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assume that this field is anchored to the dust envelope and behaves like a dipole, it should 
drop with distance cubed. However, even at ~0.5 pc from the protostar, the strength of the 
magnetic field in the jet remains comparable to that observed in the core by {30), when it 
should have decreased by more than a factor of 100 if it was merely reflecting the field in the 
ambient cloud. These results suggest that, whereas dust polarization traces the magnetic 
field associated with the ambient material close to the protostar, the synchrotron emission 
traces the magnetic field intrinsically associated with the jet. 

The polarization properties and the magnetic field configuration in the HH 80-81 jet 
are very similar to those observed in AGN jets. In AGN jets, the polarization is always 
observed either perpendicular or parallel to the axis of the jet (27). When jets are transver- 
sally resolved, the degree of linear polarization reaches a minimum towards the jet axis and 
increases towards the jet edges [e.g., Figure 4a in (31)]. Figure IB shows that the linear 
polarization in HH 80-81 is perpendicular to the jet axis, and Figure 3 shows that the de- 
gree of linear polarization increases as a function of the distance from the jet axis. In AGN 
jets, this configuration of the polarization has been interpreted as indicative of a large-scale 
helical magnetic field (27,31). 

The relevance of magnetic fields in protostellar jets has been anticipated from theoret- 
ical models of star formation [see (4,32,33) and references therein]. However, most of the 
attention has been paid to the role of the magnetic field in the launching region, and its 
relevance in the collimation and propagation at larger scales has been scarcely studied, even 
theoretically (34). Possible signatures of magnetic fields on the large scales of the YSO out- 
flows have been searched with multidimensional numerical MHD simulations of overdense, 
radiative cooling jets [see (35) and references therein]. These studies find that the effects of 
magnetic fields on the outflow are dependent on both the field geometry and the strength, 
which so far have been poorly determined from observations [e.g., (35)]. Our results for the 
HH 80-81 jet provide a measure of the strength and structure of the magnetic field in a YSO 
jet at large distances (~0.5 pc) from the driving source. Jets from YSOs produce thermal line 
emission from shocks that in turn provides direct, independent information on their physical 
conditions (density, temperature and velocity) (36) not available for jets from AGNs and 
microquasars. Our results in the HH 80-81 jet suggest that mapping of linear polarization 
in a set of protostellar jets (associated with both low- and high-mass YSOs) in combination 
with detailed theoretical modeling may lead to a deeper understanding of the overall jet 
phenomenon and of the relevance of magnetic fields in star formation. Although these radio 
jets are usually weak, these studies will be feasible with the next generation of ultrasen- 
sitive radio interferometers. Recent molecular maser observations are starting to provide 
information on the strength and direction of the magnetic field in the disks around massive 
protostars (37); observations of the dust polarization provide information on the strength 
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and direction of the magnetic field in the envelope feeding the disk in the protostellar system 
(30,38). Thus, we may be on the brink of describing the magnetic field characteristics in 
the complete envelope-disk-outflow system with a combination of observational techniques. 

Furthermore, recent theoretical works [e.g., (39)} suggest that protostellar jets could be 
a source of gamma rays. These models postulate, as a working hypothesis, the presence of 
relativistic electrons in such jets. Our detection of synchrotron emission in the HH 80-81 
jet demonstrates the presence of relativistic electrons in the jets from massive protostars, 
providing an observational ground to the theoretical conjecture and making these objects a 
potential target for future high-energy studies. 
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Fig. 1. — Images of the HH 80-81 jet region at 6 cm wavelength. (A) Image of the total continuum intensity, 
showing the whole extension of the HH 80-81 jet. The brightest radio knot (labeled as IRAS 18162-2048) 
is associated with the central protostar. The color scale [shown in panel (C)] ranges from 0.039 to 4.5 
mJy beam -1 . The rectangle marks the central region of the jet, which is shown in (B) and (C). (B) Linearly 
polarized continuum intensity image (color scale). The color scale ranges from 39 to 85 /iJy beam -1 . 
Polarization direction is shown as white bars. The total continuum intensity is also shown (contours). 
Contour levels are 40, 100, 400, 850, and 3300 /iJy beam -1 . (C) The apparent magnetic field direction 
(white bars) is superposed with the total continuum intensity image (color scale). The images were made 
by using a tapering of 20 kA in order to emphasize extended emission. The root mean square of the noise is 
0.013 mJy beam -1 , and the synthesized beam is 13" by 8" with a position angle of 2° [shown in the bottom 
left corner of panel (C)]. The (0,0) offset position corresponds to the phase center of the observations, at 
right ascension (RA, J2000) = 18 h 19 m 12.102 s and declination (DE, J2000) = -20° 47' 30.61" . 
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Supporting Online Material 

Materials and Methods 

Observations and Calibration of the Data 

Data of continuum emission at 6 cm were obtained using the Very Large Array (VLA) 
of the National Radio Astronomy Observatory (NRAO). Observations were divided into two 
runs during 2009 August, 12 and 22. Total on-source time was ~10 hours (~5 hours in 
each run). The phase center of the observations was a(J2000) = 18 h 19 m 12.102 s , 5(J2000) 
= -20° 47' 30.61". We used as phase calibrators 1751-253, in the first run, and 1924-292, 
in the second run. The amplitude calibrator was 3C 286 in both runs. Data editing and 
calibration in amplitude, phase and polarization were carried out using the Astronomical 
Image Processing System (AIPS) package of NRAO, following the standard VLA procedures. 

The polarization calibration was performed using the observations of the amplitude 
calibrator to determine the absolute polarization angle, while the observations of the phase 
calibrators were used to determine and correct the antenna-based leakage terms that produce 
instrumental polarization. The polarization of the phase calibrator will rotate in the sky 
with parallactic angle while the instrumental polarization will stay constant. We followed 
the standard VLA procedures described in chapter 4 of the AIPS cookbook. 

In order to check the quality of our polarization calibration, we measured the polariza- 
tion degree and polarization angle of the phase calibrators used in our calibration. We found 
that these parameters are consistent with those previously found by (1) and those tabulated 
in the VLA Polarization Database for the same calibrators. 

Once calibrated, data from both runs were concatenated using the DBCON task of 
AIPS. Images of all of the four I, Q, U, and V Stokes parameters were made using the task 
IMAGR of AIPS. In order to emphasize extended emission, a tapering of 20 kA was used for 
each image. The rms of the images was 0.013 mJy beam -1 , and the synthesized beam was 
13" x 8" with a position angle of 2°. 

We detected emission associated with the HH 80-81 jet in the images of the I, Q, and 
U Stokes parameters. An image of the total linear polarization emission, 
was obtained using the task COMB of AIPS. Even when the detected polarized emission is 
relatively weak (3 to 7 a), there is an element that gives it additional credibility: the polarized 
emission is observed only from the non-thermal lobes of the radio jet, with spectral indices 
of —0.8 to —0.4, and not from the much brighter central source (the thermal jet) or from 
the nearby radio sources. 
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SOM Text 

Other Possible Emission Mechanisms 

There are two other mechanisms (besides synchrotron radiation) that could produce the 
linear polarization observed in the jet of HH 80-81. However, it can be shown that they are 
too weak to be detected since their optical depths are negligible at radio wavelengths. 

A first possibility is that we are observing Thomson (i.e. electron) scattering of the radio 
emission from the central source. This mechanism would produce a polarization pattern 
similar to that observed. Assuming that part of the radio continuum emission is thermal, 
optically-thin free-free, and that the region is at a distance of 1.7 kpc (2) an upper limit of 
n e < 10 3 cm -3 is obtained for the electron density in the polarized lobes of the radio jet. 
From the observed source size of R ~ 3 x 10 17 cm (~12"), and since the Thomson cross 
section is a T = 6.65 x 1CT 25 cm 2 , we obtain an upper limit of tt = n e Ra T < 2 x 1CT 4 , for 
the Thomson opacity, clearly insufficient to produce the observed polarized emission, since 
an opacity several orders of magnitude higher is required to produce the observed polarized 
intensity in the lobes of the radio jet. A rough estimate for the radiation dispersed by this 
process is given by the flux density of the core of the jet (~ 5 mJy), multiplied by the 
fraction of solid angle subtended by the scattering region as seen from the core of the jet 
(f2/47r < 0.1) times the Thomson opacity (2 x 10" 4 ). The resulting flux density (<0.1 //Jy) 
is very small and clearly insufficient to be detectable. 

A second mechanism that could produce a polarization pattern similar to that observed 
in the jet of HH 80-81 is dust dispersion. The region surrounding IRAS 18162—2048 has 
large dust extinction. The H 2 column density is ~ 3 x 10 22 cm~ 2 (3), that translates into 
an opacity at optical wavelengths of r v ~ 30. However, the cross section for dust scattering 
scales as A -4 (4 ) an d its effect is negligible at centimeter wavelengths. 



Table SI: Flux densities and sizes of knots used to 
estimate the magnetic field strength and the energy. 





Source 


Flux Density (mJy) 


Size 


B 


E 


Knot 


Name 


20 cm 6 cm 2 cm 


(arcseconds) 


(mG) 


(erg) 


N 


15 


2.76 ± 0.14 1.62 ± 0.05 1.38 ± 0.21 


10 


0.18 


l.lxlO 44 


S 


13 


1.43 ± 0.18 0.95 ± 0.07 


10 


0.15 


7.3xl0 43 



Notes.- Flux densities and source names from (5). The 6 cm flux densities 
are consistent with those obtained in our new observations. 
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